Nonlinear perturbations for dissipative and 
interacting relativistic fluids 

David Langlois 1 ' 2 , Filippo Vernizzi 3 

1 APC (Astroparticules et Cosmologie), 
UMR 7164 (CNRS, Universite Paris 7, CEA, Observatoire de Paris) 
11 Place Marcelin Berthelot, F-75005 Paris, France; 
2 GReCO, Institut d'Astrophysique de Paris, CNRS, 
98bis Boulevard Arago, 75014 Paris, France; 
and 

3 Helsinki Institute of Physics, P.O. Box 64, 
FIN-00014 University of Helsinki - Finland 

February 5, 2008 

Abstract 

We develop a covariant formalism to study nonlinear perturbations 
of dissipative and interacting relativistic fluids. We derive nonlinear 
evolution equations for various covectors defined as linear combina- 
tions of the spatial gradients of the local number of e- folds and of some 
scalar quantities characterizing the fluid, such as the energy density 
or the particle number density. For interacting fluids we decompose 
perturbations into adiabatic and entropy components and derive their 
coupled evolution equations, recovering and extending the results ob- 
tained in the context of the linear theory. For non-dissipative and non- 
interacting fluids, these evolution equations reduce to the conservation 
equations that we have obtained in recent works. We also illustrate 
geometrically the meaning of the covectors that we have introduced. 
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1 Introduction 



Recently, we have developed PJ Oj a new formalism for nonlinear relativistic 
perturbations, based on a covariant approach [31 0] and inspired by the work 
of Ellis and Bruni 0. A key quantity in our formalism is the covector ( a , a 
linear combination of the spatial gradients of the energy density p and of a, 
the local number of e-folds, or integrated expansion along each worldline of 
the fluid. 

As we showed in PQEj, (a, for a barotropic fluid, obeys a remarkably sim- 
ple conservation equation, which is exact and fully nonlinear. In the linear 
approximation, this conservation equation reduces to the usual conservation 
law for the linear curvature perturbation on uniform energy density hyper- 
surfaces, usually denoted (. Thus, our covector ( a can be seen as a covariant 
and nonlinear generalization of the usual £■ 

It must be stressed that, although our initial motivations were related to 
cosmology, our formalism applies, in fact, to any relativistic fluid, whatever 
the underlying geometry. Thus, in the following we will not assume any 
particular type of spacetime, unless specified. 

The purpose of the present work is to clarify the geometrical meaning 
of C a , and to extend our formalism in two new directions. Whereas our 
formalism has been developed so far for a perfect fluid, we consider here its 
extension to the case of dissipative relativistic fluids. Moreover, we consider 
the possibility of having several different interacting fluids. 

The covariant formalism of Ellis and Bruni [Sj was first extended to a 
dissipative multifluid system in [HJ IZj , where both the exact and the linear 
form of perturbation equations were given for noninteracting fluids. It was 
later generalized to the case of interacting fluids in [S] for linear perturbations. 

The extension to interacting fluids is particularly useful in cosmology, 
where several types of matter coexist. As a typical example, one can mention 
the analysis of the cosmic microwave background anisotropies, which has 
been for instance considered within the covariant formalism and at linear 
order in jU] . 

In this paper we work in the framework of our nonlinear formalism 
Ej, based on covectors such as ( a . As in the single perfect fluid case, we 
show that it is possible to derive simple, exact, and fully nonlinear evolution 
equations for various covectors, constructed as linear combinations of the 
spatial gradients of scalar quantities that characterize the fluids, such as the 
energy density or the number particle, and of the local number of e-folds a. 

The equations we obtain represent nonlinear generalizations of linear per- 
turbation equations, some of which have been already studied in the context 
of linear cosmological perturbations (see, in particular, (HUE]). However, 
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our equations are covariant, i.e., independent of the coordinate system, and 
nonlinear, and it is straightforward to linearize or expand them up to second, 
or higher, order in the perturbations. Furthermore, our covariant approach 
allows to identify easily which properties belong specifically to the linear or 
second order expansion, and which properties remain valid at higher orders. 

This work is organized as follows. In Sec. El we introduce the covariant 
formalism for dissipative fluids while in Sec. H3 we discuss the geometrical 
meaning of the nonlinear covector ( a . In Sec. E]we derive an identity that can 
be used to construct conservation equations for various nonlinear covectors, 
representing the perturbations of scalar quantities of a dissipative fluid. The 
conservation equations for these quantities are derived in Sec. while in 
Sec. El we extend our formalism to interacting fluids. Finally, in Sec. we 
conclude. 

2 Covariant formalism 

In this section, we briefly review the covariant description for a dissipative 
relativistic fluid. There is a substantial literature on dissipative relativistic 
fluids, which is reviewed, for instance, in |12,.. The extension of irreversible 
thermodynamics to relativistic fluids is hampered by subtle issues. In partic- 
ular, the first extensions, due to Eckart in 1940 and to Landau and Lifshitz 
in the 1950 's suffer from non-causal behavior. An extended theory, which 
does not suffer from this problem was developed by Israel and Stewart [T3] . 
In the present work, we will not need to enter the details of these various 
formulations and, for details, we refer the reader to ^2], whose presentation 
and notation will be followed here. 

We first define the unit four-velocity of the fluid u a as the average velocity 
of the fluid particles. This means that u a is proportional to the particle 
current n a , which can thus be written as 



where n is the particle number density. 

In addition to the particle number density n, the fluid is characterized by 
local equilibrium scalars: the energy density p, the pressure p, the entropy 
S and the temperature T. In general the effective pressure deviates from 
the local equilibrium pressure so that p e g = p + U. The energy-momentum 
tensor can be written in the form 



n a = nu 



a 



u a u a = -1, 




Tab = pu a u b + (p + II) h ab + q a u b + q b u a + TC ab , 



(2) 
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where h ab is the projection tensor orthogonal to the fluid velocity u a , 

Kb = 9ab + u a u b , {h\h b c = h a c , h a b u b = 0), (3) 

and where the energy flow q a and the anisotropic stress ir ab satisfy the fol- 
lowing properties: 

q a u a = 0, n ab = n ba , ir a a = 0, n ab u b = 0. (4) 

It is useful to introduce the familiar decomposition 

V fe u a = a ab + uj ab + ^Qh ab - u a u b , (5) 

with the (symmetric) shear tensor a ab , and the (antisymmetric) vorticity 
tensor u ab ; the volume expansion, 9, is defined by 

e EE V a u\ (6) 

while u a is the acceleration, with the dot denoting the covariant derivative 
projected along u a , i.e., ' = u a V a - 

We also introduce the covariant spatial derivative, which is defined as 

D aX = h a b V bX , D aXb = h a c h b d V c Xd, etc, (7) 

where \ and \ a are a generic scalar and covector, respectively. As illustrated 
in [S], the covariant spatial derivative is particularly useful to deal with cos- 
mological perturbations in a covariant way, as an alternative to the standard 
coordinate based approach. 

The conservation of the energy-momentum tensor, 

VaH = 0, (8) 

yields, by projection along u a , the energy conservation equation, 

p + Q(p + V + n) + 7C b a V b u a + V a q a + q a u a = 0. (9) 

Using the decomposition (jSJ) and the definition of the spatial covariant deriva- 
tive 0, one can rewrite the above equation in the form 

p + 6(p + p) = V, V = - (en + 7i ab a ab + D a q a + 2q a u a ) . (10) 

The scalar quantity T> thus contains all the dissipative terms and vanishes 
for a perfect fluid. 
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In irreversible thermodynamics, the entropy is not conserved but increases 
according to the second law of thermodynamics. This can be expressed by 
the inequality 

V a S a > 0, (11) 
where S a is the entropy current. One usually writes S a in the form 

S a = Snu a + ^, (12) 

where R a is a dissipative term. As discussed in ^2] , the explicit form for R a 
varies according to the formalisms which have been introduced in the litera- 
ture. The entropy S and temperature T are the local equilibrium quantities, 
which are related via the Gibbs equation, 

TdS = d(-) +pd(-) . (13) 



This implies, in particular, 

nTS = p + Q(p + p), (14) 

where we have assumed conservation of the particle number, i.e. 

V a n a = h + Bn = 0. (15) 

Using the energy conservation equation (jlOjl . this can be rewritten as 

nTS = V. (16) 

In terms of the entropy density s = nS, this gives, using once more the 
particle conservation equation, 



s + e s = ^. (17) 



Using this equation, one gets 



In the case of a non-dissipative fluid, the right hand side is zero and the 
above relation then expresses the conservation of entropy. 
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3 Nonlinear covector 



Here we illustrate the geometrical meaning of the nonlinear covector that we 
introduced in |2] . This interpretation can easily be extended to the other 
covectors of the same form that will be defined in this paper. 

As we showed in our recent works [U |2| (see also [HJ [T5J EH] for other 
recent formulations of conserved nonlinear perturbations), a crucial quantity 
to define conserved nonlinear perturbations is the spatial gradient of the local 
number of e-folds, a, which is defined as the integration of G along the fluid 
world lines with respect to the proper time r, 

a = - f drQ. (19) 
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It follows that 

= 3d = 3u a W a a. (20) 

In [HE], we have introduced, for a perfect fluid, a linear combination of 
the spatial covariant derivatives of a and p, 

C a = D a a - -D a p. (21) 
P 

This covector is fully conserved on all scales for adiabatic perturbations, and 
can be seen as the nonlinear generalization of the usual (. In Sec. EJwe will 
rederive the conservation equation for this quantity. 

At this stage, it is instructive to give a graphical representation of ( a . 
To work with a scalar quantity rather than a covector, let us consider an 
infinitesimal vector e a which is orthogonal to the fluid four-velocity u a at 
some spacetime point q. One can then write 

c\ 

e a ( a = Aa- -Ap, (22) 
P 

with 

Act = e a V a a, Ap = e a V a p. (23) 

As shown in Fig. ^ starting from our fiducial reference point q, the in- 
finitesimal quantity Aa corresponds to the shift in a when one goes from q to 
the neighboring point q' indicated by e a . Thus q' belongs to the hyper surf ace 
Sq+ac* characterized by the constant value a + Aa. Similarly, q' belongs to 
the constant energy density hypersurface S p+ a p - Now, these two hypersur- 
faces also intersect the fluid worldline that goes through q, but in general 
the intersections differ. The quantity e a ( a quantifies, in terms of the num- 
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Figure 1: Geometric interpretation of ( a . 



ber of e-folds, the separation between these two intersection points. Indeed, 
the proper time interval between q and the intersection of the hypersurface 
S q +a« with the worldline of q is Ar a = Aa/a while the proper time interval 
between q and the intersection of the hypersurface S p +Ap with the worldline 
of q is At p = Ap/p. The difference between these two proper time intervals 
is shown in the figure and the corresponding variation of a during this time 
difference interval is e a ( a = a(Ar a — At p ). 

4 An identity for nonlinear covectors 

In this section we will derive a general identity which we will use later for 
various cases. Namely, we will show that if one starts with an equation of 
the form 



where / and g are two scalar quantities and, as before, the dot denotes the 
derivative along u a , one finds the identity 



/ + e^ = o, 



(24) 




(25) 



To show this, let us start by rewriting Equation (J24*j) as 




(26) 
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Taking the spatially projected derivative, one gets 

D a a - %Dj + -D a g = 0. (27) 
/ 9 

We now wish to invert the time derivative and the spatial gradient. In 
order to do so, it is convenient to introduce the Lie derivative along u a , C u . 
Its action on a covector is given by the expression 

£uXa = U C V c Xa + XcV a M C = U C 8 c Xa + Xcd a U C . (28) 

For a scalar, C u f = f. The Lie derivation along u a and the spatial gradient 
D a do not commute. Instead, one finds P] |2j 



D a (f) = C u (DJ) - ii a f. (29) 
Applying this identity both to a and /, one obtains 

D a a - jDJ = C u (D a a) - jC u (D a f) . (30) 
Substituting in (|27)l . one finds 

C u (^D a a - jD a fj + (j \ DJ - 3jD a g = 0, (31) 



where we have used Eq. ()24j) to rewrite the last term. Moreover, Eq. 
implies 

f a \ ' a 2 , . 

(/) (32) 

which can be used to rewrite Eq. (}3*Tj) in the form given in Eq. (|25p. 
For practical purposes, it is useful to note that, for the covectors 

= 7^ ^ = ^-7^' (33) 

which are defined as very particular linear combinations of spatially projected 
gradients, one can replace the latter by ordinary gradients and write 

CP = d a a - jd a f, = d a g - Idaf- (34) 

This is a consequence of the identity 

D a X = d a X + U a Xi (35) 

valid for any scalar quantity \- 
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5 Nonlinear (non-) conservation equations 



That one can construct conserved cosmological perturbations associated with 
each quantity whose local evolution is determined entirely by the local ex- 
pansion of the universe was shown in ^3]. There, it was pointed out that 
such a construction can be extended even in the nonlinear regime, although 
an explicit expression of the nonlinear perturbation variables was not given. 

Here we explicitly construct these variables and their evolution equa- 
tions. In particular, we use the identity f!25l) to derive conservation - or 
non-conservation - equations for various covectors which represent nonlinear 
perturbations. In this section, we discuss these equations for the nonlin- 
ear generalization of the curvature perturbation on hypersurfaces of uniform 
number density n, uniform energy density p, or uniform entropy density 
s. The corresponding covectors are obtained immediately by replacing the 
quantity / in Eq. (J25|) with n, p, or s, respectively. As in (UEj, the curvature 
perturbation is replaced here by a, the local number of e- folds of an observer 
comoving with the fluid particles. 



5.1 Particle conservation 

The particle conservation equation, 

V a (nu a ) = 0, (36) 

can be rewritten as 

n + Qn = 0, (37) 

which is exactly of the form with f = g = n. Equation (J2H|) tells us 
immediately that the particle conservation law yields, for the covector 

Cjr ] = D a a - -D a n, (38) 
n 

the equation 

£ u (^ = 0. (39) 

This result has already been given in 

This can be generalized to the case where the number of particles is not 
conserved, in which case one can write 

V a {nu a ) = -Tn, (40) 

where T is the decay rate, which is not supposed to be a constant here. 
Equivalently, one can write 

h + Qn + Tn = 0, (41) 
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which is still of the form (J2%j) with f — n and g = n + (Tn)/0. Taking the 
spatially projected gradients, one thus finds 

= — \D a V - -D a n - D a & + -D a n) . (42) 
n \ n n J 

The particle annihilation (or production) rate Tn acts here as a source for 
the evolution of the particle number density perturbation. We will discuss 
more thoroughly the case of several interacting fluids below in Sec. H3 



5.2 Energy conservation 

In the general case, if one introduces the quantity 

(3 = -| = (en + 7r ab a ab + D a q a + 2q a ii a )/Q, (43) 

then the energy conservation equation (|T0|) . becomes 

p + Q(p + p + (3) = 0, (44) 

where j3 acts as an extra pressure term, which we will call dissipative pressure. 

One can apply the identity derived in the previous section with f = p 
and g = p + p + (3. This yields directly 

C u C a = — [D aP - V -D aP + D a (3 - ^D aP \ , (45) 
P \ P P J 

in terms of the covector ( a defined in Eq. (12 lj) and introduced in |2] . 

Equation ()45|) is fully nonperturbative and valid at all scales. It holds 
for any fluid, including dissipative fluids with nonvanishing energy flow and 
anisotropic stress. It can be rewritten as 

CuCa = — (r a + s a ), (46) 
p 

where, on the right hand side, one recognizes two covectors: the nonlinear 
nonadiabatic pressure perturbation, 

T a = D a p - ^D a p, (47) 
P 

which vanishes for purely adiabatic perturbations, i.e., when the pressure p is 
solely a function of the energy density p, and a term combining the gradients 
of the dissipative pressure and of the energy density, 

E a = D a (3-^D a p, (48) 
P 
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which we will call dissipative nonadiabatic pressure perturbation. This van- 
ishes for a purely perfect fluid. 

Note that since the dissipative pressure (3 depends on the local expansion 
6, the dissipative nonadiabatic pressure perturbation E a depends implicitly 
on the local expansion. In the appendix, Sec. [XJ we discuss an alterna- 
tive formulation of the non- conservation equations which leads to evolution 
equations where the source terms do not depend on 0. 



5.3 Entropy (non-) conservation 

Defining 

P = -— = ^ (49) 

one recognizes in (fTTjl an equation of the form (ffifj) with / = s and g = s + f3. 
One can thus write the evolution equation 



£u( { a s) = — \Dj - ^D a s) , (50) 



s \ s 



in terms of the covector 

C[ s) = D a a - -D a s, (51) 
s 



which can be seen as the nonlinear generalization of the linear curvature 
perturbation on constant entropy hypersurfaces. 



6 Interacting fluids 

We now extend our nonlinear formalism to a system of interacting fluids. 
Our treatment follows the approach of PHI E], in the context of the linear 
theory, and can thus be seen as a nonlinear generalization of these works. 
In particular, we extend to the non-linear case their study of the coupled 
evolution of curvature and nonadiabatic perturbations in a multifluid system 
when energy transfer between the fluids is included. 

We work in a common "global" frame defined by a unit four-velocity 
u a . The four- velocity can be conveniently chosen depending on the physical 
problem (see [S] for a discussion on this point). In this frame, the energy- 
momentum tensor of each individual fluid can be expressed as 

= p a u a u b + p a h ab + q^u b + ql a) u a + 7$. (52) 
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In the appendix, in Sec. IB. 11 we give explicitly the transformations from the 
global frame to each individual a-fluid frame. Here, for simplicity, p a denotes 
the total effective pressure for each fluid and we will not include IP in the 
dissipative terms. 

For each fluid, the (non-) conservation of the energy-momentum tensor 
reads 

where Q^ a is the energy- momentum transfer to the a-fluid. The conserva- 
tion of the total energy-momentum tensor implies the constraint 

^Q (a)a = 0. (54) 

a 

Projecting along the four-velocity u a yields an energy conservation equation, 

p a + Q(p a +p a ) = V a + Q a , (55) 

with 

Q a = -Ql a) u\ (56) 
V a = -(7i^a ab + D a q^ a + 2q^ a u a ), (57) 

where the dot and the spatial derivative D a are now defined with respect to 
the four-velocity u a which a priori does not coincide with any fluid frame. 
If one then introduces the dissipative pressure of the a-fluid, 

P° = , (58) 

the energy conservation equation (|55|) . becomes of the form (JBJ), with 3 a — 
Q a /Q acting as an extra pressure term. This yields 

= ^(D a p a - V - a D a p a + D a (5 a -^D a p a \ 

p a V p p I 

a ( Q a \ Q a ( 6 \ 

- T D -n + w [ D - e - 7° D ° p i ■ <59) 

This is the (non-) conservation equation for 

C< o) = D a a - - a D a p^ (60) 
P 

the nonlinear generalization of the curvature perturbation on uniform density 
hypersurfaces for the fluid a. This perturbation variable associated to the 
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individual fluid a is conserved if the fluid is barotropic, P a = P a (p a ), perfect, 
V a = 0, and decoupled from the other fluids, Q a = 0. 

Note that is defined here with respect to the four- velocity u a and not 
with respect to its own four- velocity u^ a , as it was the case in ([2*T|) . This 
does not really affect the spatial gradients because they can be replaced, as 
before, by ordinary gradients. Nonetheless, the coefficient a/p a is different 
in the two cases, and one may have different definitions of Q a ) depending on 
u a . As discussed in the appendix Sec. IB.2| in the situations where the fluid 
relative velocities are small and the geometry quasi- homogeneous, as in the 
cosmological context, then the various definitions of ( a are equivalent in the 
linear theory. Their spatial projections (i.e., perpendicular to u a ) are even 
equivalent up to second order. 

To make the connection with the linear theory, in particular with [TlU lllj. 
one can rewrite Eq. (f5T)J) as 

CuC (a) = 3& + ^ _ a_ ^ )intr) + (61) 

where we have identified several individual source terms: the intrinsic nona- 
diabatic pressure perturbation, 

r W = D a p a - ^D a p a , (62) 
the dissipative nonadiabatic pressure perturbation, 

S («) = Da p" - ^D a p a } (63) 
the intrinsic nonadiabatic energy transfer, 

Q^' mtl) =D a Q a -^D a p a , (64) 
P 

and the relative nonadiabatic energy transfer, 

Si Q ' rel) = ~ [D a B - ^D aP A . (65) 

It is convenient, at this stage, to introduce the total energy density and 
pressure, 

P = E/A P = EP B » (66) 
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which are defined with respect to our unspecified global frame u a . By sum- 
ming the individual energy conservation equations (J55j) . one gets an equation 
as m (HH with 

V = Y / V a . (67) 

a 

The sum of the Q a vanishes as a consequence of the constraint ([54)1 . 

The covector ( a corresponding to p can be expressed as the weighted sum 
of the individual (( a \ 

C = Et^ (68) 

a r 

and its evolution equation is given by Eq. (j45j) . Since now the global fluid is 
made of individual fluids, it is useful to split the global nonadiabatic pressure 
T a into two terms, 

ra = r (int r)+r ( re l) ; (69) 

where the intrinsic nonadiabatic pressure perturbation is the sum of the 
individual nonadiabatic pressure perturbations, 

ri mtr) = E r i Q) - (70) 

a 

The relative adiabatic pressure perturbation can be written in the form 

where the relative entropy perturbations between the a- and /3-fluids S^ a ^ 
is defined as [10J 

s 3 ( C w _ cf) ) = _ M ^_£jjPj m (72) 
In order to derive Eq. (JoUj) from the definition (j71|) . it is convenient to use 

C^^Ca + ^E-^- (73) 
6 13 P 



Note that one can replace the spatial gradients in Eq. (JT2J) with partial deriva- 
tives. 

A similar decomposition applies to the dissipative nonadiabatic pressure 
perturbation S a , which can be written as 

S a = Si intr ) + Si rel ), (74) 
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with the intrinsic part, 



and the relative part, 



£i intr) = £si Q) , ( 75 ) 



Ea ~ ^h pp U (f) Sa ■ (76) 

Note that one could also work directly with T> a instead of /3 a , these two 
quantities being related by Eq. (|5Bj). As a consequence, one could separate 
each corresponding covector E£*) into an intrinsic and a relative part, in 
analogy with our treatment of Q^ 6 ■ 

Let us now rewrite the relative nonadiabatic energy transfer of the a-fluid 
of Eq. f)Sj) as 

e rd) = -% (D a Q - | A,p) - ^ £ £st*> , (77) 

where we have employed Eq. (J75)l for the last term. By taking the difference 
between the evolution equations (}6T| for two fluids, we finally obtain an 
evolution equation for the relative entropy perturbation, 



p q ^ y V p a p^ 3 



(78) 

The above equation represents the nonlinear generalization of the evolution 
equation for the relative entropy perturbation, established in ^T] in the con- 
text of the linear theory. 

Note that, while the total curvature perturbation ( a is sourced by the 
relative entropy perturbations between the fluids through Eqs. (fTTj) 

and dZEJ), the Ci a) 's do not appear in the evolution equation for the relative 
entropy perturbation Sff^. However, S^ a ^ is sourced by the last term of 
Eq. ([78)1 . which depends on G, i.e., on the local expansion. In the linear 
theory, this term vanishes on large scales, i.e., on scales larger than the 
Hubble radius. 



7 Conclusion 

In the present work, we have extended our covariant formulation for non- 
linear perturbations to the case of dissipative relativistic fluids, allowing for 



15 



interactions. This extension could be used to tackle more sophisticated phys- 
ical situations. Cosmology is such an example since the matter content of 
the universe is made of several fluids which can interact. 

In our approach, the important quantities are the covectors £f/\ defined 
as linear combinations of the spatial gradients of the number of e-folds a, and 
of some scalar fluid quantities /, for one or several fluids. These covectors are 
fully nonlinear and generalize the curvature perturbations on hypersurfaces of 
constant / defined in the context of linear cosmological perturbation theory. 

The non-conservation equation for ( a associated with the energy density, 
that we obtain when the fluid is unperfect, can be derived and written in a 
form analogous to that found in our initial formalism with a single perfect 
fluid, provided that the pressure is modified such as to include dissipative 
effects. 

For several interacting fluids, we can also define fully nonlinear quan- 
tities that generalize analogous quantities introduced in the linear theory. 
Remarkably, as in our initial formalism with a single perfect fluid, the evolu- 
tion equations that govern these quantities "mimic" those found in the linear 
theory, with the advantage that they are covariant, i.e., independent of any 
coordinate system. Thus, it is straightforward to linearize them or expand 
them at second, or higher, order in the perturbations. 

Moreover, our fully nonlinear approach allows to identify which property 
is specific to linear (or second) order and which property will remain valid 
at all orders. For example, in the multifluid case, there are a priori several 
nonlinear generalizations of the curvature perturbation, which depend on the 
choice of the reference frame. In the cosmological context, they all reduce to 
the same quantity at linear order. 

A Alternative formulation of the non-conser- 
vation equations 

In order to apply our general identity (J25J) to the case of a dissipative fluid, 
we have introduced the dissipative pressure j3 which depends on the local 
expansion 0. In doing this, we have obtained equations like (|45jl where 
the source covectors, on the right hand side, depend implicitly on the local 
expansion. As we show below, there exists an alternative formulation, which 
leads to evolution equations where the source terms do not depend on 0. 

Instead of writing the non-conservation equations in the form ()24|) . an 
alternative possibility is to start from an equation of the form 

f + Qg = h. (79) 
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Taking the spatially projected derivative of the above equation and re- 
placing by 3d, one gets 

DJ + 3gD a a + 3aD a g = D a h. (80) 

For the first two terms one can invert the gradient and time derivative by 
using the identity (J3HJ). One obtains 



Cu (D a f) + SgC u (D a a) = D a h - 3aD a g + hu a , (81) 

where we have used Eq. ()79|) to simplify the terms proportional to ii a . Using 
again Eq. (|7U|) to replace 3d on the right hand side, this relation can be 
rewritten after some straightforward manipulations in the form 

C u [D a a + -^jDjJ =-^5 {9°-f ~ f D -9 + hD "9 ~ gD a h) + ^-u a . (82) 

We have thus derived an evolution equation for the covector 

((J)=D a a+^DJ, (83) 

which is a linear combination of the spatial gradients of a and /, but which 
differs from the covector defined in f!33j) . when h in ()79j1 does not van- 
ish. These two quantities have also different geometrical interpretations. In 
particular, unless h = 0, does not generalize the curvature perturbation 
on the hypersurface of uniform /, as does. However, one advantage of 
this alternative formulation is that the number of e-folds a, or its proper 
time derivative, the expansion 0, appears only on the left hand side of the 
non-conservation equations. On the right hand side, one only finds quantities 
characterizing the fluid, although the acceleration now appears. 

This new identity can be applied to the case where the number of particles 
is not conserved, in which case one can write 

n + Qn = E, (84) 



where H is the production rate of particles. Eq. (|82|) with f = g = n and 
h — H then yields directly 

^ = "I (^ " ¥ " <) • < 85 > 

The new identity also applies to the non-conservation of energy in the 
dissipative case, with f = p, g = p + p and h = V. This yields directly, for 
the covector 

C ( a p) = D a a + ^-pjDaP, (86) 



17 



the evolution equation 

£u( { a p) = - 3 , * \pD aP - pD a p + VD a {p + p) 

-(p + p)D a V} + ^—jU a . (87) 

The first two terms on the right hand side are proportional to the covector 
r a defined in flH|). 



B Multifluid system 

B.l Frame dependency of the energy- momentum ten- 
sor 

In a multifluid system, the energy-momentum tensor of each fluid a can be 
written in the form 

# = P a u^ + p«h% + q^ + q ( b a) u^ + (88) 

where the unit four- velocity of the fluid a - e.g., the rest frame of 

the particles of the fluid a - and is the projector on the hypersurface 
orthogonal to u^ a , 

Quantities with the tilde are written in the a-fluid frame. 

Each fluid can have a different frame u^ a . It is preferable, however, to 
work in an arbitrary global frame defined by a common unit four- velocity u a . 
This is related to the four-velocity of each individual fluid by the relation 

u { : ] = 7a (u a + «<«>) , u a v^ = 0, (90) 

with the relativistic factor 7 a = (1 — g ab v^v^)^ 1 / 2 and v^ a a peculiar ve- 
locity. The energy momentum tensor of each individual fluid can be rewritten 
in the global frame as 

= p a u a u b + p a h ab + q^u b + qi a) u a + ir%> . (91) 

The quantities in the (untilted) global frame can be expressed in terms of the 
tilted quantities and of the peculiar velocities. The corresponding expressions 
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are given by [T7] and read 

P a = 7«p Q + (7«-l)p Q + 2 7Q gt ) ^ )a + ^ ) ^ ) V a ) fe , (92) 

q { a ] = ll(p a +PnV^+laqi a) +lai a) V^ b (vi a ~>-U a ) 

+7t^v^ b - nffv^ b v^ e u a , (94) 

= iI^ + p^^^^-Ik^v^) 

+7. (i°M a) + ^ - ^ a) « w j 

+jr^ ) v (a)c v^ d u a u b - -h ab f}^v^ c v^ d . (95) 

3 

These equations allow to reexpress in the global frame defined by u a , the 
energy- momentum tensor of each fluid a written in their rest frame u^ a . 

Note that at linear order in the perturbations, if the fluid relative veloc- 
ities are small and the geometry quasi-homogeneous, as in the cosmological 
context, then the energy density, pressure and anisotropic stress of the fluids 
do not depend on the frame. Only the energy flow depends on the peculiar 
velocity v^ a . 



B.2 Frame dependency of 



The covector Q a ) defined in Eq. (IfiOj) for a fluid a depends on the choice 
of the global four- velocity u a . However, for each individual fluid, one can 
define a nonlinear generalization of the curvature perturbation on uniform 
density hypersurfaces, only with respect to its own four-velocity i/ a K This 
is defined as 

go.) = D (a) a _ ?&!^ D (a) (96) 

where 

D^ X = h^ b V bX , (97) 

and is defined in Eq. (JHHJl . Note that one can replace the spatial gradients 
with standard partial derivatives in Eq. ([96)1 . 

Let us now assume the situation where the relative velocities between the 
fluids are small and let us consider a common reference velocity u a so that 
the differences 

e (a)a = u (a)a _ (gg) 
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are small. One can then treat the e^ a as perturbations and the expansion 
of the coefficient in Eq. yields, at first order in e^ Q ) a , 

(99) 



_ .(a ay n a, i 
" Va " ~ _ + - e ( a ) a ( 

Inserting this result into the definition of Q a \ one gets 



= d a a - ^^daP - C { a a) + U a ^\i a \ (100) 

Consequently, the projections of the two covectors and Q Q ) perpendicular 
to u a are equivalent up to first order in e^ a . 

One can go one step further when the spacetime is approximately homoge- 
neous, as in cosmology, because the spatial gradients can then be considered 
as perturbations as well. Then, the difference between and Cjf^ * s °f 
second order. This implies that spatial parts of the two covectors, i.e., their 
projections perpendicular to u a , are equivalent up to second order. 
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